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Section 1 

SUMMARY. CONaUSIONS m RECOMMENDATIONS 
1 .1 BACKGROUND 

The objective of this test program was to evaluate the performance of a star 
tracker which uses a G.E. ST-256 CID as the focal plane detector. The CID is 
an array of 256 x 256 pixels which are 20 x 20un in dimension. The tracker 
used for test was a breadboard tracker system developed by BASD. Unique ac- 
quisition and tracking algorithms are employed to enhance performance. A pat- 
tern recognition process is used to test for proper image spread function ard 
to avoid false acquisition on noise. A very linear, high gain, interpixel 
transfer function is derived for interpolating star position. 

The Nikon F-2.8 lens used in the tracker has an EFL of 100 mm. The tracker 
has an FOV of 2.93 degrees resulting in a pixel angular subtense of 41.253 arc 
sec in each axis. It should be noted that the CID used in this tracker is a 

research device procured from GE. The device was furnished with peripheral 
circuits packaged in breadboard fasfvion. It ivas essentially the first of its 
kind and had not been tested beyond a level necessary to insure it was opera- 
tional when received by BASD. The tracker's projected performance was based 
entirely on data taken from a previous 128 x 128 CID fabricated and tested by 
G.E. 



The test program of this report, conducted with funds provided under MSFC con- 
tract No. NAS8-34263, was designed to evaluate only the tracker system per- 
formance. No detailed CID characterization tests were conducted under this 
program. 


1.2 APPROACH AND EXTENT OF TEST 

Appendix A is a. copy of the test procedure (BASD No. 2332-101) used for the 
program. This procedure presented a star to the tracker in a circular pattern 
of positions; the pattern was formed by projecting ^ simulated star through a 
rotatable deviation wedge. 
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Two diameters were used for the circle* a small diameter of approximately 4 
pixels and a large diameter of approximately 14.5 pixels. 

Five small circle locations were tested with the CIO temperatures at 0*C, 10*C 
and 20*C and with star magnitudes of approximately 1.0, 2.8 and 6.0 Mv. Thus, 
nine small circles of data were taken at each of the five locations. Each da- 
ta circle consists of 60 data points taken at approximately 6 degree intervals 
of wedge rotation. Each data point represents the mean of 5 readings. 

Each small data circle sampled approximately 65 different pixels and used 16 
different interpixel transfer function cycles in each axis. 

The large circles were taken only at 0*C and 2.8 Mv. The number of data 
points and readings were the same as those used for the small circle read- 
ings. Each large circle sampled approximately 720 different pixels and uses 
30 different transfer function cycles in each axis. 

Table 1-1 shows the CIO extent for pointing accuracy. 

The pointing accuracy test was automated to take 60 data points for a circle 
in less than 10 minutes. The raw pixel data for the track cycle was stored on 
floppy discs. The computer then determined positioir and error for the data 
relative to a best fit circle. All data was normalized to the pixel dimension 
of 2Q-.jm. 

Further tests determined readout noise. Noise Equivalent Displacement (NED) 
during track, and spatial noise during acquisition by taking related data and 
reducing it. The standard deviation for a number of readings at a fixed posi- 
tion was used to determine readout noise and NED. Coarse and fine maps were 
developed using the tracker acquisition mode. These maps were examined to de- 
termine the raaximun spatial noise as it, applies to the coarse Md fine acquis- 
ition modes. 



Table 1-1 

EXTENT OF ACCURACY DATA 



Small 

Large 



Circles 

Circles 

Total 

No. of locations for circle 

5 

5 

10 

No. of circles of data 

45 

5 

50 

No. of different pixels used 

325 

3,600 

3,92_5 

No. of different transfer function cycles 

160 

300 

460 

No. of position measurements 

2,700 

300 

3,000 

No. of readings taken 

13,500 

1,500 

15,000 
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Extensive exploratory tests and analyses identified potential sources of error 
and methods of improving performance. 


The test results are sunmarized in paragraph 1.3 below. Further detail on the 
test approach and results are provided in Sections 2, 3 and 4 of this report. 
Volume II contains a complete set of the data that was taken. 


1.3 SUMMARY OF TEST RESULTS 


All position, dimensional and accuracy data has been normalized to the 20u 
pixel dimension and can be related to angle by the scale factor of 41.253 arc 
secs per pixel. 

Accuracy data related to an arbitrary reference on the CIO is all inclusive. 
No attempt has been made to separate error sources such as instability, test 
equipment, quantization, image uncertainty, etc. 

Tracking data was taken using an LED source with a center wavelength of 656um 
and the tracker lens set at F-2.8. 

1.3.1 Interpixel Transfer Function Linearity 

The linearity of the interpixel transfer function was evaluated over 'several 
cycles and results indicated that there is less than + q. 005 pixels-RMS devia- 
tion from a straight line. This performance is equivalent to that predicted 
by BASD’s analyses. 

1.3.2 Interpixel Signal Scale Factor 

The scale factor (Signal vs. Displacement) for pixel-to-pixel transfer is ap- 
proximately 1.7 times the total image signal per pixel displacement. The the- 
oretical limit is 2.(1. The expected scale factor was approximately 1.9 times 
the image signal per pixel displacement. 
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The scale factor is inversely proportional to image diameter. For this syten 
the minimum “apparent image diameter" appears to be defined by the pixel-to- 
pixel cross talk in the CID. 

1.3.3 Readout noise 


A random readout noise of 278 e" per pixel was observed for the CIO focal 
plane and is consistent with G.E.'s prediction and is reduced to approximately 
35e" per update by the multiple read Non Destructive Read Out (NDRO) process 
developed by G.E. 

If the system bandwidth is equal to the pixel data rate of 5000 Hz, the read- 
out noise constant is 3.93 e" per nh 

1.3.4- Spatial Noise For Coarse Acquisition 

The double read subtraction process for coarse acquisition takes approximately 
4 seconds miminum and results in an integration time of 2 seconds which can be 
increased by adding integration time between reads. A maximum spatial noise 
of 10^e“ was observed. 

BASDiS-desrign is based on a minimum star signal of 3.3 x lO^e" per second. By 
adding two seconds integration to the acquisition time this signal should rise 
above the noise level and result in an acquisition time of 6 seconds for a 
minimun star signal level. 

1.3.5 Spatial Noise For Fine Acquisition 

Using an integration time of 0.5 second fine acquisition maps revealed a worst 
case spatial noise- of lO^e". This is below the minimun star signal level and 
fine acquisition can take place in less than 1 second without false attempts. 
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Acquisition Performance 


The acquisition modes were not optimized beyond the point required to acquire 
stars for accuracy testing. Extensive testing in this area was considered to 
be inefficient at this time; however, acquisition was readily achieved for the 
1 Mv to 7 Mv star range used fn this test program under the basic conditions 
of paragraphs 1.3.4 and 1.3.5. A significant number of false acquisition at- 
tempts were observed for stars dimmer than 6 Mv, especially at the 20*C CID 
operating temperature. 

The pattern recognition process effectively rejected acquisition on spatial 
noise when thresholds were below noise level. 


1 .3 J Position Accuracy Performance 

For stars brighter than 3 Mv the accuracy performance was + 0.0178 pixels 
RMS. There was no significant variation in accuracy using CID temperatures 
ranging from 0*C to 20*C. The errors for large circles were slightly larger 
than for small circles, a degradation attributed to changing focus resulting 
from nonorthogonality of the CIO plane to optical axis alignment. 


The accuracy for the 6 -Mv star roagattude did vary with CID temperature chan- 
ges. The results from all the 6 Mv star data are as follows: 


CIO TEMPERATURE 

0*C 

10*C 

20* C 


ACCURACY 

+ 0.0282 RMS 

t 0.0266 RMS 
+ 0.0465 RMS 


These results are inconsistent with BASD's previous analyses which showed 
tracker accuracy would be limited by pixel-to-pixel varizrtions in dark cur- 
rent. Therefore, all data should show degraded accuracy with increasing CID 
temperatures. The errors for the 6 Mv star magnitude do change with tempera- 
ture but do not follow the predicted trend and are higher then expected. 
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1.3.8 Sources Of Error 

The analysis indicated that accuracies of +.005 pixels to +.015 pixels RMS 
should be achievable for CIO temperatures of 0*C to 20*C respectively. 

Review of data for stars brighter than 3 Mv does not reveal a trend of de- 
graded accuracy with increased CIO dark current. The transfer function line- 
arity error (Subsection 1.3.1) meets the criteria for the projected accuracy. 

Another error source, not previously identified, was suspected and an exten- 
sive evaluation and.exploratory test program was conducted by BASD to identify 
the error sources. 

The conclusion is that the accuracy for bright stars could be limited by pre- 
sent mechanisation protLlems rather than by uncertainties in the CIO (except 
for bad areas in the CID). A solution to this problem has been defined but 
not yet implemented and tested. The problem and solution are briefly summar- 
ized below. 

• The CID design provides access to a 4 x 4 block of pixels with one 
address. 

• BASD's mechanization requires sampling a dark pixel for use as a 
reference for common mode noise rejection. 

• In order to conserve time, the dark pixel is sampled from the same 
4x4 block used to sample the star signal pixels. This results in 
sampling “dark" pixels adjacent to the pixels on which a star is 
imaged. 

• Sufficient pixel-to-pixel cross talk exists in the CID; therefore, 
a pixel adjacent" to the star image also receives star sigriaT. 
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• The cross talk into the dark reference pixel results in biasing the 
transfer function for varying star positions which could result in 
position bias errors. 

The problem can be minimized by advancing the data block when taking dart da- 
ta, and thereby assuring that the dark pixel is never sampled adjacent to the 
signal pixels. Although this takes more tin>', it is estimated that the re- 
sultant increase in NED will be less than lOX. 

Other potential sources of error that have not been investigated are: 

• Pixel-to-pixel response variations 

• Pixel point spread response function and variations 

• Effects of the injection pulse 

• Test equipment errors (non-flat wedges) 

• CID misalignment (non-orthogonal alignment) 

• Image distortion 

• Chromatic effects of Cid and lens 

1.3.9 Noise Equivalent Displacement (NED) 

The-NED-was tested~u s i ng a: 5^9~Hv star with varying update time. Results are 
tabulated below in Table 1-2. 

The NED performance closely follows the predicted pattern in that for snort 
update times it is dominated by readout noise, and for longer update times it 
is dominated by signal and dark current shot noise. 
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Tdole 1-2 

NED VS. UPDATE TIME (Tu) 
Mv = +5.9 


NED (RMS Pixels) 


Tu (seconds) 

x 

y 

.080 

.0695 

.0431 

.150 

.0215 

.0208 

.250 

.0086 

.0083 

.500 

.0043 

.0037 

.800 

.0028 

.0025 

1.300 

.0032 

.0029 

Tracking Rate 




Initial attenipts to conduct the accuracy test resulted in frequent loss of 
track while moving the wedge, and consequently the star image to a new posi- 
tion. Due to this problem and an excessive long acquisition time for dim 
stars, star pattern recognition logic and a portion of the adaptive rate con- 
trol were implemented. Thereafter, few loss of track events were experienced, 
and when they were, the star was reacquired with little time lost in taking a 
circle or data. 

The adaptive rate control allows for tracking without degradation in accuracy 
to rates of 0.5 pixels per update. There is an associated lag in position da- 
ta and a slight decrease in the transfer function slope; however, both are 
predictable. 

The tracker should maintain track to rates equivalent to 1 pixel per update. 

The test was not set-up for precision measurements at high rates and the 
tracking accuracy could not be verified; but an open loop test was conducted 
to verify track was maintained at 1 pixel per update. This was accomplished 
by driving the star in a large circle pattern at maximum wedge rate and esti- 
mating the pixel rate of the image. The update time was increased for an 
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equivalent rate of 1 pixel per update. The rate used was approximately .5 
pixels/second. 

1.3.11 General Performance of the ST-256 CID 


Detailed characterization tests for the CID were not part of this program; 
however in the coarse of conducting the tracker test some general observations 
were made. 

• CID environment? and operating time 

The ST-256 was exposed to uncontroTTed handling and environments typical 
of mo?t R & D programs. During the development and troubleshooting pro- 
cess the focal plane was removed and reinstalled in the optical head many 
times. The CID sealed glass coverplate was inadvertently cracked during 
a microscope inspection early in 1981 and was subsequently removed. A 
dry N2 purge was added tc provide an inert environment for the exposed 
detector. Unobserved loss of the Nj supply during a 0*C test resulted in 
frost formation on the CID chip. 

The unit was transported between facilities at BASD. It was shipped to 
Washington D.C., May 1981 and subsequently subjected to low level radia- 
tion to evaluate noise susceptability. 

BASD has operated the CID on an almost continual basis since its receipt 
in July 1980. Estimating its on-time at SOX, approximately 1700 hours 
have accumulated on this device without noticeable performance degrada- 
tion. 

• CID Noise 

The star tracker readout noise is slightly less than the 300e“ predicted 
by G.E. 
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Spatial noise was not evaluated, but, it appears to approximate G.E.'s 
i predictions. Refinements to the acquisition and track algorithms will 
^ provide more insight into average spatial noise limitations. 

' Some abnormal local dark current generation was observed as was expected. 

However, the frequency of occurrence should be low relative to the total 
CID area. 

This CID appears to have several high current areas or "hot spots". 
These hot spots were noted in false acquisition attempts and, in some 
cases, in false tracking when thresholds were set low. No quantitative 
measurements were made except for the area where accuracy tests may have 
been performed. 

This particular device also exhibits a large electrically coupled noise 
pulse immediately following injection. The anomally was observed by G.E. 
prior to delivery and was attributed to an error in the chip mask. The 
amplitude of this pulse can be reduced by a mask change. Common mode 
noise rejection tn the tracker design minimized the effect on tracker 
performance. The most noticeable effect was ^ reduction in the usable 

dynamic range of the A/D converter. The same problem relates to the 
fixed pattern noise evident in the CID. 

• CID dark current 

The CID dark current was continually monitored in accordance with the 
procedures of Appendix A. The measurement made using the dark reference 
pixels errs somewhat since the pixels contain some signal (subsection 
1.3.8). 

The measured mean dark current is approximately ID^e", 2.4 x lO^e" and 5 
X lO^e" per pixel -secon<Lat jO*C, lO'C and 20'C respectively. These le- 
vels were obtained from dim star tracking and are consistent with the ex- 
pected results. 
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• CID Responsivity 

Star magnitude is determined by using the track signal pixels and comput- 
ed in accordance with Appendix A procedures. The tracker Mv computation 
is made using the spectral response curves obtained from a G.E. report on 
the 128 X 128 CID. A test determined that the computed Mv magnitude and 
that measured with a photometer agree within i-5 Mv for stars in the 
range of 1.0 Mv to 6.0 Mv. This is close enougn to lend credibility to 
the predicted CID response. 

1.4 CONCLUSION 

The trackers performed with an accuracy of approximately i *02 pixels RMS at 
0*C. The tracker's scaling is 41.253 arc sec per pixel resulting in an accur- 
acy of t 0*82 arc sec over the 2.93 degree FOV. 

The accuracy of + ,02 pixels, or 1 part in 12,800 for the total FOV, is im- 
pressive when compared with sensors using past technologies. However, it does 
not meet the degree of accuracy BASD feels the tracker can ultimately achieve. 

From the results of this test and previous studies, it is concluded that the 
CID star tracker remains a primary candidate for continued development and 
performance-evaluation. A potential error resulting from the existing .mechan- 
ization as discussed in subsection 1.3.8 can be corrected. Other potential 
sources of error remain to be investigated. Since our projected linearity 
(better than +0.005 pixel) for the interpixel transfer function was verified 
and since accuracy is not yet limited, at least for brighter stars, by the CID 
dark current, it appears that further improvements can and should be made. 

1.5 RECOWiENDATIONS FOR FUTURE WORK 

A variety^ of areas require further study, development and evaluation. The ex- 
tent of effort associated with each depends on the objectives to be accom- 
plished. The effort suggested below is directed at defining the CIO tracker 
capabilities and shortcomings in order to eliminate technical risk from future 
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flight development programs. The reconmended future efforts related to both 
the CID and the tracker system design are summarized below. 

1.5.1 Proof of System Concepts 

BASO has completed extensive analyses and studies related to: 

• Acquisition 

• Automatic Gain Control (AGC) 

• Adaptive Rate Control 

• Multiple Star Tracking 

• Tracking Accuracy 

The completion of the following recoirmendations will result in a higher degree 
of confidence relative to the expected performance in each area. 

Improved Accuracy 

A potential systematic source of error has been identified in this test 
program and can be eliminated with a software change. It is recommended 
that the change be implemented in the existing equipment and that the 
test be repeated to verify maximum tracking accuracy for a tracker with- 
out compensation or calibration. Other potential sources of error should 
also be investigated under this task. 

Improved Acquisition and Rate Tracking 

A portion of the acquisition logic and adaptive rate control were mecha- 
nized to complete this test program. Additional functions must be added 
to optimize performance in both areas. It is recommended that these 
functions be added and evaluated in the existing system. 
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Multiple Star Tracking and AGC 

Thresholds and gains can be readily selected for proper acquisition and 
track of a single star. However, for multiple star tracking the signal 
from the related sources may be several orders of magnitude apart. To 
acquire, select, and track these sources in a reasonable amount of time, 
it is necessary to automate the acquisition ranging and gain of the track 
loop. BASD has studied this problem and formulated various concepts to 
achieve the desired performance. It is reconmended that these concepts 
be refined and implemented- in the existing hardware for evaluation. 

1.5.2 CIO Study and Eva^luation 

Recoimiendations for study in the CID performance and the readout process ap- 
plicable to tracker acquisition and track modes-are as follows. 

CID Characteristics and Specification 

Several CID's should be fully characterized and evaluated in order to es- 
tablish confidence in their performance in future trackers. A specifica- 
tion for the CID oerformance parameters should be written reflecting the 
results of this evaluation program. 

CID Design and Readout Approach 

The CID design and the readout approach should be thoroughly reviewed 
with the following objectives: 

• Elimination of injection pulse 

• Reduced pattern noise 

• Reduced readout noise 

• Improved data and processing rates 

• System Simplification 

Any significant improvements should be breadboarded and tested. 
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— Mecham'cal/ThermaT Design 

General packaging concepts should be studied In order to establish a 
stable design which can withstand launch and space environments. The de- 
s^lgn should provide adequate cooling for the detector and a means for 
heat transfer of cooler power. A mechanical model should be fabricated 
and environmentally tested to ensure these requirements are met. 

1.5.3 General Study Items 


• Optics 

An analysis of the full spectrum of CID spectral sensitivity for typical 
lens performance and CID point response should be conducted to relate 
tracker performance versus star color. Whenever possible, cromatic veri- 
fication tests should be conducted using breadboard hardware. 

• Parts Status and Power Estimates 

A detailed study should be conducted on parts for use in a flight unit 
and should include plans to qualify or change nonstanxlard parts. A de- 
tailed power estimate should be made based on a typical design for flight 
application using reasonably acceptable parts. 
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Section 2 

TEST APPROACH AND EQUIPMENT 
2.1 GENERAL APPROACH 

A compound reflective/refractive system which projects a light source to in- 
finity was used for the star sinulation. The star light was projected through 
a motorized, compound, rotatable deviation wedge assembly. As the wedges are 
rotated in the drive assembly, the star appears to move in a circular path. 
This circular locus of the star image is the positional input to the star 
tracker optical assembly for position accuracy tests. 

The star tracker optics consist of the ST-256 CID detector mounted at the fo- 
cal plane of a Nikon F-2.8 camera lens. A Cromemco Z-80 microprocessor devel- 
opment system is interfaced to the star tracker optics which performs the ac- 
quisition and track logic functions. The star simulator, wedge drive asseobly 
and star tracker optics are mounted on an isolated test surface plate to oro- 
vide maximum mechanical stability. 

A block diagram of the test setup is shown in Figure 2-1 



Figure 2-1 Block Diagram of Overall Hardware System 
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2.2 STAR SIMULATOR 



The star simulator projects a light source from an LED to infinity. Monochro- 
matic light was used to eliminate optical problems associated with ordinary 
white light when using a coimercial lens. The LED output is a green light of 
approximately 656nX wavelength. 

The intensity of the simulated star was varied by an external power supply. 
Three star intensities used gave a range of about five star magnitudes. The 
EFL of the simulator is 61 cm. An aperture diameter of 0.051 mm was used in 
front of the LED which resulted in an angular subtense of 17.25 arc sec. (less 
than half the pixel dimension). A photograph of the star simulator is shown 
in Figure 2-2. 

2.3 DEVIATION WEDGE DRIVE 


The wedge drive assembly is located in the optical path between the star siai- 
lator and star tracker optics. The assembly contains two independent sets of 
counter rotating gears. Either set may be independently driven by a D.C. mo- 
tor coupled to the gears. Four deviation wedges with adaptors are provided 
with the assembly. The deviation angle of each is approximately 12 arc min. 
Counter Tot aling tbe^ weitges produces a linear deviation of “the star light as a 
function of wedge rotation. A single rotating wedge simulates circular motion 
of the star. 


Figure 2-3 illustrates the concept. 
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Figure 2-3a. Two Counter Rotating Wedges Produce Linear Motion 



Figure 2-3b. Single Rotating Wedge Produces a Circular Motion^ 
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When two wedges are attached to one gear, indexing their deviation vectors 
produces a desired resultant deviation. Both wedges rotating together result 
in a circular motion with an adjustable angular diameter of 0 to 24 arc min. 


A scribe line on each wedge indicates the direction light will be deviated. 
Each wedge housing has 132 teeth, on its periphery at intervals of 2.727 de- 
grees. These teeth, along with the wedge scribe, allow for course indexing of 
wedge pairs to achieve a desired resultant deviation. A lock for the gears 
farthest from the motor is provided to select the set- which is driven. The 
two sets are coupled through planetary friction drive gears. 


A one percent continuous turn servo potentiometer is coupled to one gear of 
either counter rotating set. The coupling gear ratio is 1:5 so that with one 
complete gear rotation the potentiometer rotates through five cycles. Each 
pot rotation corresponds to 72 degrees of gear and associated wedge rotation. 

The control box which drives the wedges can be set to either incremental steps 
of wedge rotation (manual or computer initiated) or free^ continuous run (set 
manually). The speed of the rotation is manually variable from zero to a nax- 
imum. The wedge rate control is open loop in either case. Figure 2-4 is a 
photograph of the wedge drive assembly. 

2.4 STAR TRACKER 


The CID is a 256 x 256 pixel detector array with a pixel size of 20 x 20 ]ga. 
The array is designed for parallel row readout and includes integral shift re- 
gisters, The- CID signals are processed using double read nondestructive read- 
out (NDRO). The process reduces fixed pattern and readout noise. 

The effective focal length of the Nikon lens is 100 mm giving each pixel an 
- angular subtense of approximately 41.253 arc sec. The array, therefore, pro- 
vides a square field of view of approximately 2.933 degrees on an edge. The 
CIO imager is mounted on a board which contain? the hardware circuitry for 
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Figure 2-4 Wedge Drive Assembly 
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reading out the array. The board also has the CID cooler circuit on it. Tne 
analog CID signals are directly digitized by an A/D converter to usable form 
for the microprocessor. 

Figure 2-5 is a photograph of the CID and circuit board. Figure 2-6 is a 
photograph of the complete tracker optics. 

2.5 COMPUTER HARDWARE 

The CID tracker data is acquired and processed in software. Complete control 
and analysis of this system is performed by the Z-80 microprocessor through a 
Crotnemco development system. The interface box consists of several computer 
hardware cards with interfaces to the CID, terminal, line printer, and plot- 
ter. The terminal displays real time data as the star tracker acquires and 

tracks stars. Data is stored on floppy disks and then printed out by the line 
printer. The plotter then plots the acquired data. 

Figures 2-7 and 2-8 are photographs of the computer hardware and peripheral 
test equipment, 

2.6 GENERAL TEST CONDITIONS 

During tests, the room lights were off and other light sources were mini- 
mized. Operation was conducted at ambient room temperature and a hood was 
placed over the tracker, wedge assembly, and simulator to provide a mere 
stable temperature environment by minimizing convection currents. The CID had 
a constant N2 purge of about 20 lb. while the system was running and about 5 
lb. when the tracker was off. The data was taken at three different CID tem- 
peratures. Temperature variation is accomplished through use of the T-E cool- 
er upon which the CID is mounted. 
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Figure 2-5. CID Board and Lens 
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Figure 2-7 Complete Hardware 
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Section 3 
TEST RESULTS 


3.1 NOISE TEsrr 

3.1.1 General Considerations 

Noise is represented by basic random or spontaneous fluctuations resulting 
from the physics of the devices and material which comprise the electrical 
system. The performance of the CID tracker depends upon the noise generated 
within the system; which is categorized as temporal noise (random noise} and 
pattern noise (spatial noise) are of particular concern. 

3.1.2 Tempora.1 Noise 

Temporal noise is a totally random signal and the two major sources are John- 
son noise and shot noise. Johnson noise is the dominant temporal noise source 
at high rates and is proportional to the square root of the bandwidth. Shot 
noise is proportional to the square root of the signal, and in this case, the 
signal can be light or dark signal /current. Shot noise typically becomes dom- 
inant at high temperatures or long update times. 

3.1.3 Pattern Noise 


The two types of pattern noise which are of concern are "fixed" pattern noise 
which repeats with each reading and "variable" pattern noise which is not re- 
peatable. The latter may vary according to time, temperature, and signal 
level. 

The major sources of fixed pattern noise in CID image sensors are transistor 
switching interference, array photolithographic variations end bias charge 
variations. The variable types of pattern noise result from spatial varia- 
-tions^ in dark -current and-responsivity. 
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The device used in this test program had another noise source resulting from 
an error in the CIO mask design^ This noise occurred as a pulse coupled into 
the signal line immediately after injection and had a long decay time and gave 
the CIO the appearance of having a shaded spatial time dependent response. 
This problem can be eliminated or minimized in future devices by a correction 
to the mask design. 

3.1.4 Noise Suppression Techniques 

The temporal noise is reduced by repeated CID nondestructive readings or up- 
dates. This process effectively reduces the Johnson noise by a factor equiva- 
lent to the square root of the number of readings taken and is used in^Jxith 
the acquisition and track modes of operation. 

Fixed electrical pattern noise is minimized by the double read, nondestructive 
method of readout. Two consecutive sets of read samples are differenced to 
cancel the repetitive spatial noise interference, ideally leaving the signal 
behind. 

Dark current pattern noise can be reduced by cooling the CID to minimize the 
mean dark current. The variations (spatial noise) in dark current and respon- 
sivity are inherently low in the CID. Further reduction is accomplished by 
averaging the effects over several pixels in the readout-process^ 

The effect of the injection pulse was minimized using common mode rejection 
techniques. A pattern recognition process implemented in the tracker results 
in overriding acquisition on noise. The noise primarily effects the time for 
acquisition and tracking accuracy. 

3.1.5 Noise Tests 


Tests were conducted to define the temporal and pattern noise evident in the 
CID tracker. 
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3. 1.5.1 Temporal Noise Data 

The temporal noise test was conducted on a block of 16 pixels. 128 nonde- 
structive readings (NDRO's) were taken for each pixel. The second 64 readings 
were subtracted from the first 64 to eliminate fixed pattern noise. This pro- 
cess was repeated a hundred times for each pixel in the 4x4 pixel block. 
Thus, the (spatial) fixed pattern noise is minimized and random or temporal 
noise is left. The test was repeated at O'C and 10*C. The mean and standard 
deviation for the sets of a hundred readings are shown in Table 3-1 and 3-2. 
- The mean of the readings decreases in value with time as they are read from 
pixel #1 to pixel #16 because of the injection pulse. 

The average standard deviation at 10*C is 17.59 bits. The scale factor is 
179. 5e” per bit, so that the temporal noise for 128 readings is 3157. 4e“ 
RMS. The average standard deviation for the 0* data is 3134e~ RMS. Taking 
the average of the 10* and 0* data, there are 3146e" RMS for ^he 128 read- 
ings. Since the noise increases as the square root of the number of readings, 

3146 278el 

the noise due to one reading is: or ^adout • 



3-3 
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Table 3-1 

TEWORAL DATA FOR 1& PIXELS AT 0*C 


Pixel # 

Average of 100 Readings 
X (of Signal) 

S (Standard Deviation) 

1 

58a.7475 

20.0954 

2 

584.4343 

25.2458 

3 

588.6263 

16.3386 

4 

594.3939 

20^900 

5 

185.6162 

14.0391 

6 

185.7980 

16.4248 

7 

187.8283 

16.549J3 

8 

183.3131 

14.5059 

9 

117.4242 

16.4112 

10 

120.8485 

18.4515 

11 

114.7172 

18.6449 

12 

123.1010 

16.2195 

13 

84.3939 

16.7087 

14 

77.6465 

16.3715 

15 

83.9091 

17.4228 

16 

85.0404 

15.2569 
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Table 3-2 

TEW>ORAL DATA FOR 16 PIXaS AT 10* C 


Averaqe of 100 Readings 

Pixel # X (or Signal) S (Standard Deviation) 


1 

700.22 

17.3523 

2 

7D.06 

16.8612 

3 

697.1 

16.0407 

4 

709.05 

17.2060 

5 

239.38 

16.4984 

6 

302.71 

22.4423 

7 

286.05 

18.0903 

8 

289.18 

18.3361 

9 

214.19 

16.7680 

10 

223.02 

18.3809 

11 

220.60 

16.8810 

12 

222.18 

16.7244 

13 

190.20 

15.6870 

14 

189.61 

16.5015 

15 

197.51 

22.1088 

16 

189.75 

15.4111 
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278 

The noise constant normalized to frequency is: kn = 

kn = 3.93e~/ /Hz 

This assumes the bandwidth is equivalent to our sample rate of 5000 Hz. 

3. 1.5. 2 Pattern Noise Data 

Volume II of this report, (Section 1.1) of the test data contains coarse and 
fine maps. Data is presented in hexadecimal units. Each reading on the 
coarse map is the sum of a block of 4 x 4 pixels. Each reading on the fine 
map is that of a single pixil. The scaling for these maps is 179.5e"/bit. 
The magnitude of the readings is the conibined dark current and system noise. 
The effect of the injection pulse is evident on the coarse map. A graph show- 
ing the magnitude of the coarse map reading against time is plotted in Figure 
3-1, which shows the effect of the injection pulse on the data. Tne corre- 
sponding data is presented in Table 3-3. 

The parameter of significance in the coarse and fine maps is the difference in 
signal for adjacent readings. By comparing adjacent readings star acquisition 
is determined. Some significant observations are summarized in Table 3-4. 

3.1.6 Conclusions For Noise Test 

The observed random readout noise was slightly less than the 300e"/sample re- 
ported by G.E. The pattern noise had not been previously tested for the me- 
thod implemented in this tracker. 




Table 3-3 

COARSE MAP READINGS 
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Set # Block # t(msec) Pulse (e") xlO^ Set # Block # t(msec) Pulse (e") xl9^ 


1 

.41 

7.399 

4 

16 

85.66 

2.870 

2 

.82 

6.361 


32 

92.22 

2.949 

3 

1.23 

6.254 


48 

98.78 

3.086 

4 

1.65 

6.277 


64 

105.34 

2.601 

3 

2.06 

6.256 


16 

111.9 

2.710 

6 

2.47 

6.223 


32 

118.5 

2.815 

7 

2.88 

. 6.038 


48 

125.0 

2.980 

8 

3.30 

5.830 


64 

131.6 

2.520 

9 

3.71 

5.521 

6 

16 

138.2 

2.705 

10 

4.12 

5.441 


32 

144.7 

2.750 

12 

4.94 

5.302 


48 

151.3 

2.981 

14 

5.77 

5.179 


64 

157.8 

2.468 

16 

6.59 

5.119 

7 

16 

164.4 

2.675 

20 

8.24 

5.272 


32 

170.9 

2.707 

24 

9.89 

5.338 


48 

177.5 

2.985 

28 

11.54 

5.369 


64 

184.0 

2.430 

32 

13.18 

5.087 

8 

16 

190.6 

2.646 

36 

14.83 

4.386 


32 

197.2 

2.646 

40 

16.31 

4.933 


48 

203.7 

3.003 

44 

17.80 

4.936 


64 

210.3 

2.513 

48 

19.70 

4.979 

9 

16 

216.8 

2.574 

52 

21.42 

5.001 


32 

223.4 

2.600 

56 

23.07 

4.805 


48 

230.0 

2.955 

60 

24.72 

4.681 


64 

236.5 

2.457 

64 

26.40 

4.141 

10 

16 

243.1 

2.571 

16 

32.96 

3.542 


32 

249.6 

2.556 

32 

39.55 

3-,606 


48 

256.2 

2.996 

48 

46.14 

3.570 


64 

26.8 

2.470 

64 

52.73 

2.933 

16 

64 

422.0 

2.357 

16 

59.30 

3.060 

32 

64 

-844.0 

2J95 

32 

65.90 

3.107 

18 

64 

1265.0 

2.233 

48 

72.50 

3.226 

64 

64 

1687.0 

2.276 

64 

79.10 

2.725 






tr {time to read one block) = .41 msec for this coarse map data 


Map scaling = coarse reading (decimal) xl79.5e” 

Fvt 

Results shown are the difference of two samples taken approximately 2 seconds 
apart. 
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Fiqure 3-1. Injection Pulse vs. Readout Time 
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Table 3-4 

SIGNAL DATA ON COARSE AND FINE MAPS 
Fixed Pattern Noise Data Summary: 



MAP 

CID 

TEMP 

TR (TOTAL 
READ TIME) 

Ti (TIME 

BETWEEN 

READINGS) 

SIGNAL A 

SIGNAL (MAX) 

Coarse (Dark) 

0‘C 

3.85 sec 

0.125 sec 

218272. Oe '/block 101417. 5e’ 

Fine (Dark) 






Location 






20 

20 

O’C 

32 msec 

460 msec 

4846.5e'/pix 

3410. Oe" 

15 

15 

0*C 

32 msec 

460 msec 

5923.5e"/pix 

5564. Oe- 

10 

10 

0*C 

32 msec 

460 msec 

6103.0e'/pix 

7439. Oe' 

20 

20 

16*C 

32 msec 

460 msec 

18309. Oe-/pix 

5385.0e- 

15 

15 

16*C 

32 msec 

460 msec 

19924. Oe'/pix 

5205. 5e~ 

10 

10 

16X 

32 msec 

460 msec 

20104. Oe"/pix 

9693. Oe- 


The maximum differential signal, A signal, between adjacent blocks is what 
sets the lower acquisition threshold limit. A signal should ideally be zero 
since fixed pattern noise and dark current have been cancelled. The differen- 
tial signals shown above are the worst occurrence of residual pattern noise, 
dark current variation and Johnson noise in that location of the CID array for 
the condition shown above. 

The worst case differential signal (or pattern noise) for adjacent blocks in 
the coarse map at 0*C is approximately 10 electrons. Setting the coarse ac- 
quisition threshold above this level results in no false acquisition attempts. 
Our system analysis has considered a signal of approximately 10^e"/sec for the 
minimum star. A worst case geometry of the star location relative to a block 
could reduce this to 3 x 10^e~/sec (or 1/3 of the minimum star). 
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The minimum coarse acquisitloa time is 4 seconds and yields an effective inte— 
gration time of 2 seconds. Therefore, approximately 2 seconds of integration 
time must be added to acquire the minimum star. The total acquisition time 
would be 6 seconds. 


The fine acquisition noise is more favorable. Note that at 16*C, location 10, 
10 has a maximum differential signal (noise) of lO^e” for approximately a half 
second integration time. Using the v*orst case star signal of 3 x 10 e /sec, a 
fine acquisition time of less than one second should result in no false acqui- 
sition attempts. The total acquisition time for both coarse and fine maps for 
a minimum star should be less than 7 seconds. 


3.2 ACQUISITION PERFORMANCE 


The acquisition mode for the tracker has not been optimized and no extensive 
tests were conducted. However, with one set of conditions (fixed thresholds, 
NORO's and integration time) the stars over five magnitudes of intensity were 
readily acquired in less than 10 seconds of time. The pattern recognition 
concept implemented in the tracker effectively rejected areas of high pattern 
noise which had previously caused problems in acquiring the dimmer stars. Now 
stars can be acquired when the star signal is significantly less than the pat- 
tern noise. 


3.3 POINTING ACCURACY TESTS 

3.3.1 General Procedure 


The pointing accuracy test was performed in accordance with BASD procedure # 
2332-101 (see Appendix A). The test verified the position interpolation algo- 
rithm developed by BASO. The design goal was to achieve an error of less than 
112 of a pixel deviation for the outpu*: star position from the input position 
of the star image. Tests were conducted over small areas involving greater 
then 50 pixels and large areas involving more then 700 pixels. 
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The fine or small field accuracy tests were performed by rotating the star im- 
age in a circular path of fixed radius. The small circles were centered at 
five different locations on the CID array- Position data was recorded for 
three different star intensities and at three varying CID temperatures at each 
location. 60 positions around the circles were sampled five times each re- 
sulting in 300 data points per small circle. The mean of these five samples 
is the position output for each of the 60 positions. 

The large field pointing tests were used to determine the relative pointing 
accuracy over a larger portion of the CID. The star image was rotated in a 
circular path with a radius of about 14 pixels. Position data was collected 
for five large circles centered at five different locations. See Figure 3-2. 



LARGE FIELD ACCURACY SMALL FIELD ACCURACY 

Figure 3-2. Pattern For Test Data 

The circular motion of the star image was obtained by the wedge drive assembly 
for both the small and large field accuracy tests. Depending upon wfiich wedge 
was used, either large or small circles could be traced. (See test equipment 
and procedure Section 2-0 for more details.) The voltage input to the star 
simulator was varied to obtain the three different light intensities. The CID 
temperature was controlled by a potentiometer setting on the CIO electronic 
board which is connected to the computer interface box. 

The data was collected and processed by the Z-80 microprocessor systea. The 
four constants needed fo’' the position interpolation mechanism were computed 
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iteratively in software. The constants were computed for each new location to 
avoid variances due to focus changes. The interpixel transfer function was 
assumed to be a straight line (Output vs. Displacement) over each half pixel 
of displacement. The four constants are used to normalize the output to the 
pixel dimension. The data was analyzed by applying the least square method to 
fit the 60 position points to a perfect circle. Appendix B discusses the 
arithematic used. 

The center location and radius was computed for the circle. The displacement 
of each data point from this circle is the interpolation error for that posi- 
tion, The standard deviation of all data points is the interpolation error 
for that circle of data. These errors were computed twice, the first compu- 
tation using all the position data accumulated and the second using only the 
position data deemed valid. Invalid data may have resulted from spurious 
noise spikes while reading or computing position and were identified by a high 
standard deviation of the five samples taken for that position. Individual 
data points consisting of five readings that had a standard deviation greater 
than one percent of a pixel were discarded in the second computation. 

With the help of software, the equipment design allowed taking 300 points in 
about 10 minutes. Transcribing the data onto the floppy disk was the single 
most time consuming factor. This automation reduced the possibility of test 
errors from room temperature variation and mechanical instability of the sys- 
tem. 


3.3.2 Pointing Accuracy Test Data 

The pointing accuracy test data is summarized in Tables 3-5 thru 3-7. The raw 
data for a CID temperature of 0*C is contained in Section 4.0 of this vol- 
ume. The total set of data is included in Volume II of this report. Note 
that all position,^ dimensional and error data was normalized to the pixel di- 
mension of 20um. The scaling for this system is 41.253 arc sec/pixel. 
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Table 3-5 
SUMMARY DATA FOR 
POINTING ACCURACY TESTS 

ORIGINAL PAGE 
OF POOR QUAUTY 
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Table 3-6 
SUMMARY DATA FOR 
POINTING ACCURACY TESTS 

ORIGINAL PAGE IS 
OF POOR QUALITY 
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Table 3-7 
SUMMARY DATA FOR 
POINTING ACCURACY TESTS 

ORIGINAL PAGE IS 
OF POOR QUALITY 
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3. 3. 2.1 Description of Sunmary Sheet (Tables 3-5 through 3-7). 


1) Disk/File #: 


All data is stored on floppy disks and identified 
by the file number shown. 


2) CID-TEMP: The temperature of the CID is set by a potention^e- 

ter on the CID electronic board. Temperature shown 
is in degrees centigrade. 


3) Focus REF DIM: The lens focus is monitored by a micrometer for 

each new circle position. Reference focus dimen- 
sion is shown in mils. 


4) Star magnitude: This is computed by an equation developed in Sec- 

tion 5.32 of BASD procedure No. 2332-101 (see Ap- 
pendix A). 


M. = la > ! <>’ WRQ (Ir - t Z I U) , j.si 

S {f/n)‘ 

4.465 X 10^ is a constant for the interpolation Algorithm 
NDRO = total number of samples in first and second read 
Tr = total time spent taking the first and second read 
Ti = integration time between the first and second read 
S = signal in e" (179.5e’/bit) 
f/n = lens f-number setting. 


Tr = tr NDRO 


tr * time for one read of track-pattern 
(12 pixels) = 2.9 X 10“^ sec 


5) Effective Charge Integration time (Tieff) = .5 Tr + Ti: This is the 

effective charge integratioa time for a. g^Lven_total 
read time and integration time (Ti) between both 
readings (See Figure 3-3). 
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6) Update Time (Tu) 


is the total time taken to obtain a reading for 
tracking information. (See Figure 3-3). 



Tieff 


Figure 3-3 Update Time 


The number of readings taken in the read time Tr 
during the track cycles. 

The dark current generation due to the effective 
charge integration time (Tieff) 


NORO (Tr + 2 Ti) sec-pix 

where 

NDRO total # of readings for first and second reads 
S IS the output signal in e” 

(Invalid data is excluded from the following column) 


7) Track NDRO: 


8) Id (e'/sec-pixel): 


as 
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9) Center Location (x,y): The center coordinates of each of the small or 

large ctrcles. 

10) Radius (pixels); The mean radius of the circle computed from the 

least square method of fitting the sixty points to 
a circle. 

11) Interpolation Constants: Computed iteratively in software to 

(Kx, Kx^, Ky, Ky') normalize the transfer functions to the pixel di- 
mensions. 

12) Standard Deviation: The deviation of position points from the least 

squares fit circle positions. 

13) Points with Error >: A histogram of all sixty points for the number of 

points with deviations in the following ranges, .05 
to .04, .04 to .03, .03 to .02, .02 to .01, and .01 
to .00. 

3. 3. 2. 2 Description of Computer Data Sheet (Test Data Section 4.0 of This 
Volume and Volume LI) 

Tach data printout contains information on one circle TJf 300 ijosition points. 
Typical data includes the number of NDRO's, integration times and thresholds 
for coarse, fine and track modes. Also, each time a position is sampled, the 
corresponding temperature of the CIO, wedge angle, dark current generation, 
star magnitude and positions are recorded. 

Computer Printout Details: 

A computer sheet containing the following details is printed for each circle. 

• NDRO's is the total number of samples in the first read only, for 
coarse, fine and track. (The same number is used for the second 
read but is not printed.) 
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• Integration Delay is the time, in seconds, between the first and 
second read for coarse and fine acquisition and track. 


• Thresholds is the minimum threshold, in electrons, for acquisition 
in coarse, fine and track modes. 


• Constants are the four constants (Kx, Kx', Ky and Ky') used to nor- 
malize the transfer function to the pixel dimension. 

Tabulated Data 


• Temperature is the CID temperature given in degrees centigrade. 

• Angle is the deviation wedge angle potentiometer outp'it. The wedge 
angle rotates 6* for each position sample. The potentiometer to 
wedge gear ratio is 5:1, hence for one rotation of the continuous 
pot, the wedge rotates through 72*. 

• Dark Current is the dark current 3 enerated in the effective charge 
integration time (Tieff). Units are in e”/sec-pix. 

• Star Magnitude is computed using a formula developed in Section 
5.32 of procedure No. 2332-101 (Appendix A.) See Paragraph 3. 3. 2.1 
item 8. 

• Reference defines the coordinates for the origin from which the in- 
terpolated star position will be computed. It also defines the lo- 
cation of the star to the nearest half pixel. 

• Position defines the coordinates of the computed star position. 
These are the sum of the reference location and the interpolated 
position of the stari -Each data point is the mean of 5 pos-it-ion - 

readings. 
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Standard deviation refers to the variation in the five readings de- 
fining the position. 


Reduced Data 


• Delta Ri are the errors for the 60 data points shown in the first 
tabulation. Reading the data from left to right corresponds to the 
sequence of the coordinate data tabulated above. The error shown 
is the straight line deviation of the coordinate position from a 
best fit circle. 

• XO, YO, RO are the center coordinates and radius for the best fit 
circle as determined by the least square method. 

• Standard deviation is the standard deviation of the sixty data 
points from the best fit circle. 

• Histogram Deviation is a tabulation of the number of Delta Ri that 
fall within the five ranges; .00 to .01, .01 to .02, .02 to .03, 
.03 to .04, .04 to .05. The reduced data is computed a second time 
with bad data removed as described in paragraph 3.3.1* 

3. 3. 2.3 Plots 

The plots of small and large circles consist of all 60 Oata points. Points 
23, 27 and 29 did not plot because of a software problem. The standard devia- 
tion printed below the plot is computed from valid data points only. The fol- 
lowing are also included: the disk/file number of the corresponding data, 

lens f/n number, focus reference dimension, the four interpolation constants, 
star intensity and the temperature of CID during test. 

Although the plots are not useful for quantitative evaluation of performance, 
they were extremely useful in troubleshooting the hardware/software and inde- 
terraimng constants. See test data. 
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3-3.3 Star Magnitude Calibration 

A star magnitude calibration was conducted to verify the star magnitude equa- 
tion 

( 4.465 X 10^ KDRO ^.51 

s (f/n)^ 

The star magnitude varies as a logrithmic function so that 

(2.512)”'' (1 X 10 -^ 2 ) = watts/cm2 

where here a OMv star is estimated at 1 x 10"^2 ^tts/cm^: 

A UDT photometer measured the flux output of the star light for an intensity 
at a zero magnitude star. A factor of 10 corresponds to approximately a fac- 
tor of 2.5 change in ({i|y|y. Since 10 percent of maximum star was measured to be 
4.92 Mv, the maximum star corresponds to 4.92 - 2.5 ~ 2.4 Mv and oiniraum star 
(1 percent of maximum) would correspond to 7.4 Mv. These star magnitudes 
closely approximate the magnitude (to about +.5 Mv) computed by the software 
formula for Mv above. 

3.3.4 Noise Equivalent Displacement (NED) and Update time (Tu) 

The noise equivalent displacement (NED) is the random error in star image cen- 
troid position resulting from detector and circuit noise. This evaluation 
measured NED as a function of update time (Tu). The noise equivalent ais- 
placement analysis is based on (RMS) noise which decreases as the number of 
samples read. The number of samples are given by: 

NDRO = yJ - 


wnere 

NDRO IS the number of samples for the first and second readout 
Tu is the update time. See Figure 3-3. 
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Ti is the integration time between reads 

tr is the read time per sample of the track subarray (2.9 x 10"^ sec). 


The tes^t involved recording approximately 25 samples of star position for a 
stationary star in the track mode. Data was collected for the conditions of 
zero degree centigrade (CIO temperature), minimum star magnitude (~ +5.8 to -^6 
Hv) and minimum integration delay time. Different number of NDRO's were used 
for each set of 25 samples to vary the update time as snown by the equation 
above. 


The results of tnis test are shown in Table 3-8. 


3.3.5 Hap for CID Region with High Current Generation 

Some of the raw data plots show highly erratic data points as a result of ab- 
normally high local current generation in the CID. This is particularly ev:- 
dent in the small circle data taken in the first position. 

A fine map was taken in this region (see Volume II of this report) in an at- 
tempt to quantify the anotnmaly. The data shows differential dark current in 
the pixels of approximately 4 x 10 e and 10 e at 10*C and 20*C respective- 

5 

ly. A-typical map -would show variations of less than 10 e- under the sene 
conditions. 

3.3.6 Transfer Function Linearity 

Analysis has shown that BASD's interpixel transfer function non-linearity 
should not exceed +i percent for image diameters greater than one pixel. 
Several sets of data were analyzed and plotted against a straight line func- 
tion. Figure 3-4 shows the typical deviation from a straight line for this 
data and shows a maximum deviation of approximately t ^.5 percent. 
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Figure 3-4. Typical Deviation From a Straight Line For Data 
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Table 3-8 
NED Data 


Ti = = ~ .4 nisec 
Hv = +5.9 


NED 


NDRO 

Tu (sec) 

X 



y 

10 

.080 

.0695 



.0431 

ID 

.150 

.0215 



.0208 

30 

.250 

.0086 



.0083 

4E 

.400 

.0062 



.0045 

60 

.500 

.0043 



.0040 

75 

.600 

.0043 



.0037 

SA 

.700 

.0028 



.0046 

9C 

.800 

.0028 



.0025 

AD 

.900 

.0030 



.0031 

□ 

1.000 

.0032 



.0029 

FF 

1.300 

.0032 



.0029 

A p^lot of^ HED vs. 

update time- 

is shown 

in 

Figure 3-5. This plot shows that a 

noise equivalent 

displacement 

of less 

than 1 

percent of a pixiel can be met 

with update times 

as low as .23 sec. 

In 

this 

case, the best performance was 

achieved with Ti minimum and Tu of .8 sec. 

for 

a NED of 0.25S of a pixel. NED 

for the X axis is 

higher than 

NED for 

the 

y axis because of a misalignment in 


this particular tracker. Proper alignment should yield a uniform NED for both 
axes. 
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Evaluation Of Error Sources 


An extensive evaluation and exploratory test program was conducted by BASD in 
order to identify the error sources. The conclusion drawn from this program 
fs that the accuracy is limited by our current mechanization rather than by 
uncertainties in the CID (except for bad areas in the CIO). A solution has 
been defined, but not yet Implemented and tested. This problem and solution 
are briefly summarized below. 

a) The CID design provides access to a 4 x 4 block of pixels with one ad- 
dress. 

b) BASD's mechanization requires sampling a dark pixel for use as a refer- 
ence for common mode noise rejection. 

c) The dark pixel is sampled from the same 4x4 block which is used to sam- 
ple the star signal pixels to conserve time. In some cases this results 

in sampling "dark" pixels adjacent to pixels on which the star is imaged. 

d) Sufficient pixel -to-pixel cross talk exists, internal to the CIO, so that 

a pixel adjacent to the star image also receives star signal. 

e) The cross telk into the dart reference pixel results in biasing the 
transfer function differently for varying star position and results in 
position biaj errors. 


The problem can be minimized by advancing the data block when taking dark data 
so that the dark pixel is never sampled adjacent to the signal pixels. Al- 
though this requires more time, it is estimated that the resultant increase in 
NED will be less than 10 percent for a given update time. 

Other potential sources of error remain to be investigated. 


i 

i 
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3.3.8 Conclusion for Pointing Accuracy Test 

The accuracy for stars brighter than 3 Hv is approximately +0.0178 pixels over 
a CIO temperature range of 0*C to 20*C. The accuracy for all data at 0*C is 
better than +0.02 pixels. The accuracy for the 6 Hv star degrades with CID 
temperature. 

From the data of Tables 3-5 through 3-7 the cross talk is quite evident. Note 
that the dark current is highest for the bright star in all sets of data. 
This -results from the higher amount of signal coupled with the dark reference 
pixel from which dark current is computed. The dark current variation for the 
raw data sets is also indicative of this effect. 

The temperature dependency of the 6 Mv star data is much greater than expected 
(if It is totally due to CID dark current variations). This would indicate a 
mean pixel -to-pixel variation in dark current of approximately 5 percent. 
This is excessive compared to the 2.2 percent measured on the 128 x 128 CIO. 
It is not known at this time whether or not the cross talk problem may be am- 
plifying these effects. Before further conclusions can be drawn, further in- 
vestigation needs to be maoe into this problem and into the problem of high 
current generation spots on the CID. 

The-cros^s- talk- was previously taken into account on a qualitative basis in de- 
veloping algorithms for a linear transfer function. However, it was not con- 
sidered when the cormon mode noise rejection was introduced. 

Except for the cross-talk bias effects, the CID tracker performed as expected 
within the limits to which the acquisition and track algorithms were refined 
for this test program. The test approach and associated hardware -developed by 
BASO proved efficient and accurate. There is no indication that test errors 
are influential in the accuracy test data. However, further investigation in 
this area remains to be completed. The test software developed under this 
contract results in a hard copy of all pertinent data. Except for the. time 
required to compute constants, a circle of 60 data points can be taken and an- 
alyzed in approximately 20 minutes with an estimated accuracy of approximately 
0.1 seconds of arc. 
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Section 4 

POINTING ACCURACY DATA AT 0*C 


4.1 0*C TEST DATA 

The data in this volume has been limited to the accuracy test at 0*C because 
it is of primary interest and the total set of data would not conveniently fit 
in this section. 

Volume II contains a full set of data. The printout is described in Paragraph 
3.3.2.2. 
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Section 1 
INTRODUCTION 

This document defines the procedures, that will be used to conduct the test 
program described in BASO's proposal No. P078 titled, "Test Program for Bread- 
board Stellar Tracker System," 

The test that will be conducted are: 

1. Inter-pixel Transfer Function 

■ 2. Small Field Pointing Accuracy 

3. Large Field Pointing Accuracy 

4. Rate-Tracking Performance 

5. Acquisition Sensitivity 

6. Tracking Sensitivity 

7. Tolerance to CIO Defects 

The accuracy tests of one through four and-±est seven will be conducted using 
a system of light deviating wedges to control the angle of a light beam emina- 
ting from a star simulator. 

The wedge drive system allows linear deflection by counter rotating two wedges 
-or conical deflection by rotating a single confound wedge. The angular change 
that is introduced will be determined by monitoring wedge rotation angle via a 
potentiometer geared to the wpdge. The observed angular change will be 
determined by the star tracker output. Differences between the two define the 
star-tracker accuracy. 

The procedures described in the following sect' .'s are designed tc provide a 
test accuracy of better than +o.23s or less than • 5 percent of the CID. pixel 
dimension. Further discu:sion_jf test errors is presented in Section 5. 

Star intensify changes required for Tests 5 w.J_6 will be accomplished through 
the use of neutral density filters. 
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Section 2 

TEST EQUIPHEFIT AND SET-UP 

The general test equipment antJ set up Is shown in Figure 2-1. The salient 
features of the equipment Is described below. 

2.1 STAR SIMULATOR/AUTOCOLLIMATOR 

The star simulator is a compound reflective/ refractive system that projects a 
light source to infinity. The Effective Focal Length (EFL) of the optics is 
approximately 48 inches. 

The source is an illuminated aperture approximately 0..Q023 inches (0.00584 cm) 
-in ■dtam'^er that is illuninated by a tmgsten filament lamp. 

The angular substense of the projected source is therefore approximately 
lOS. 

The spectral content of the simulated- star wflT not be defined for this pro- 
gram. However, it will be maintained constant- for all tests by operating the 
lamp at a constant current. Int^rsity will be varied by introduction of neu- 
-tral-dentity-fi Iters. 

The star simulator-may also he used-as an-autocul limator through the use of an 
eye piece in the optical path. The FOV in the autocollimator mode is ±2. 5m. 
This feature will be uied for alignment of the wedge rotation axis and for 
coarse alignment of the tracker. 

2.2 DEVIATION HEDGE DRIVE ASmiiBLY 

The drive assemoly contains two independent sets of counter rotating gears. 
Either, set may be independently driven- bjt i. DC- motor that -is coupled -to— the 
gears. 
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Four deyiation wedges with adaptors are provided with the assembly. The 
deviation angle of each Is approximately 12m. Th& wedges mey be adapted and 
oriented In several configurations for varying effects. Some of the alterna- 
tives are as follow: 

a. Attach a wedge to each of the gears. Each counter rotating set 
will the yield linear deviation as a function of wedge rotation. 

This configuration can be used to provide Independant horizontal 
(x) and vertical (y) displacement of the star In the focal plane of 
the tracker. 

b. Attach one wedge to one gear so that circular motion is provided In 
the focal plane.. 

c. Attach two wedges to one gear and index their deviation vectors to 
achieve a desired resultant deviation. Rotation now results In 
circular motion with an adjustable-angular diameter. (0 to 24m) 

d. Attach two wedges to each gear of one set. Linear, or ellptical, 
deviation of adjustable amplltude-occurs when counter rotated de- 
pendlng-on the indexing of each wedge pair. 

Each wedge- has a scribe line on its edge Indicating the direction light will 
be deviated. 

Each wedge housing has 132 teeth on its perlfery at Intervals of 2.727 de- 
grees. These teeth along with the wedge scribe, will allow- for coarse In- 
dexing of wedge pairs to achieve a desired resultant deviation. 

A lock for the gears farthest from the motor Is provided to select the set 
--that is dr iv en ^ When disengaged that set will-drive-. When-engaged the- other 
set win drive. The two sets are coupled through planetary friction drive 
gears. 
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A one percent linear potentiometer is provided idtir the system and may be 
coupled to one gear of either counter rotating set. The coupling gear ratio 
Is li5 so that, the potentiometer rotates five times for one rotation of the 
gear. Each pot rotation is therefore 72 degrees of _gear rotation. This 
yields an equivalent accuracy in reading gear position of +0.2 percent or 
0.72 degrees, neglecting backlash and tooth variations. The test process de- 
fined in paragraph 3 and 4 reduces the effects of this error to less than O.IS 
01^0.0025 pixels^as^seen by the tracker. 


2.3 STAR TRACKER 


The star tracker lens 1s a Nikon f/2.0, 50 irm camera lens. 

The ST-256 CIO is mounted at the focal plane of the lens. The CIO is a 
256 X 256 pixel detector array with a pixel size of 20 x 20 wn. 

The ER of the lens is 100 mn which results in a pixel angular subtense of ap- 
proximately Al'.TSSS. The array therefore provides a square field of view of 
approximately 2.9335 degrees on an edge. 

The CIO tracker data is acquired and processed through the use of a Comemco 
Z-SOTutcroprocessor development set. 

A has-ic acquis-lttoa ancUtrack. mode is mechanized so that the FOV is searched 
until a "star" \s found. The "star" is then tested to insure that it is not 
just a faulty pixel. If it is a faulty pixel, acquisition search is 
continued. If it is not, it will be tracked until it leaves the FOV or is 
lost by other means such as too high a rate. When it is lost the acquisition 
mode is reinitiated. 

BASO's tracking algorithm does not require position information so that star 
positfon is Ttot computed in the basic track loop. Star position can however 
be determined from the signals generated while tracking. 
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Ilvis- Is^ -accomplished by recording "the tracking data and computings position 
after-the-fact» 

Critical ^tracker parameters such as integration time and number of NDRO's are 
controllable from the terminal. It is also possible to obtain individual 
pixel data for detailed evaluation of CIO defects. 

2.4 RAW TEST DATA 

The flexibility of the star tracker development equipment allows taking and 
prlnting-auvariety of_data.^ 

The minimum raw data that will be furnished for this test is tabulated 
below. Other data wfTl be furnished as deemed necessary during the test pro- 
gram. 


Test Data 

a. Test set-up data, hand written (each new set-up). Typical informa- 
tion-will -include: 


• General wedge-orientations and configuration 

• Star simulator current 

• K7D^ -fftters-being used 

• Star tracker temperature 

• Lens: f/n 

• Test objectives 

• Number of NDRO's used 

• Integration time 

• Update-time 

-w Time and date 

b. Tracker data, printed (each data point) 




X-axis output 
Y-axis output 


2-5 



TR81-04 



• Star Intensity output 

• Relative wedge (gear) angle output 

• CIO temperature 

The tracker data will be printed In decimal and will have, attached to It, the 
related test set-up data. 

Data analysis to determine performance will- be conpleted In accordance with 
Instructions described In the procedures of Section 4. All position data will 
be handled In terms of the pixel linear dimension rather than Its angular sub- 
tense. Conversion to angular position can be readily made as desired. 
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Section 3 
TEST PROCEDURES 

The procedures described In this section describe only the processes required 
to obtain the desired test data^ Parameters such as star Intensity or lamp 
current, CIO temperature^ lens f/n, number or NDRtTs,. integration time, up- 
dated time, pot voltage, etc. shall be optimized by the operator unless speci- 
fic requiremnts are calTedrTor In the procedures, 

"This Information along with other pertinent information about test conditions 
shall be recorded on each set of printed data either by writing It on the data 
sheet or by attaching a separate data sheet, 

The procedures herein are established withr-the -intent to minimize test equip- 
ment and process errors. Should It become evident that alternate procedures 
or data would enhance the test program, the procedures will be altered. 

All data associated with 1 nterpoTatin^ position shalT be taken with e single 
compound wedge and the rotation shall be confined to one direction to minimize 
backlash errors. 

3.1 INTER PIXa TRANSFER FUNCTION, SMALL FIELD POINTING ACCURAa, 

EFFEaS OF CIO PARAMETERS AND RATE TRACKING ACCURACY 

The objectives of this test are to: 

a. Define the constants required for true position Interpolation 

b. Evaluate relative accuracy over relatively small areas of the CIO 
(several pixels) 
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c. Evaluate effects of CIO dark current as related to position ac- 
curacy 

d. Evaluate rate tracking accuracy 

The set-up and data taking for these test is combined because similar position 
data is required to evaluate all associated parameters. - • 

This test evaluates relative performance over several small portions of the 
total CIO FOV, each of which is several pixels in diameter. 

3.1.1 Set-up Requirements 

a. Assemble two wedges to gear of the front set and the potentiometer 
and two wedges to one gear of the aft set. 

Adjust the front set so that the deviation vectors (scribes) ere 
coincident. 

Adjust the aft set so that the deviation vectors (scribes) are 
-separated-by approx-imately-65- teetlv-on-the-housings. 

Identify- the front set -as "coarse" and-the aft set as “fine."- 

This configuration will result in a circular scan of the star image for either 
wedge set when rotated. For the fine set the diueter will be approximately 
3.5- pixels and for the-coarse, approximately 70 pixels. 

b. Adjust the integration time between readings (NDRO's) to- zero or 
the minimum that can be achieved. 

c. Adjust the star brightness for the maximum possible while maintain- 
ing proper tracking performance. 
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d. Align the star sioulator so the star is near the center of the 
tracker FOV (typically +50 pixels) 

e. Place a plane parallel mirror against- tht-wedge drive gear and ad- 
just the drive assembly until the return image ir in the star simu- 
lator/autocollimator FOV. 

f. Adjust the CIO temperature to approximatley 0*C. 

3.1.2 Data Requirements 

a. Obtain tracker position wedge angle and intensity data for 60 posi- 
ttons~of“^he fine wedge- -uniformly spaced at intervals of six 
degrees wedge rotation. The tolerance for wedge postton settings, 
as read from the pot, shall be +2.5 degrees and all settings shall 
be acquired by rotation in the same direction^ 

Identify data as “small field pointing, max star, 0*C“ (note,, if 
temperature is not 0“C, record actual temperature)-. 

-b^ -Decrease the star intensity by 50 percent -and increase the integra- 
tion time until ouput intensity is the same as that obtained in (a) 
above. Use approximatley the- same wedge rotation angles. Repeat 
the test or 3J..2 (a) above and identify data as “small field 
pointing, 50 percent star, 0*C“. (Record integration time.) 

c. Oecrease-the star another 50 percent and- increase the integration 
time for the same intensity read out as above. Identify data as 
“small field pointing, 25 percent star, 0*C.“ (Record integration 
time.) 

d. Restablish the conditons of 3.1.2 (a) except for the CIDS tempera- 


ture. 
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Repeat the test of 3.1.2 a, b, and c for the CIO temperature set at 
10*C and then again for the CIO temperature set at room tempera- 
ture. 

e. Obtain the pixel outpouts from a 12x12 subarray centered on the 
locus of positions used for the tests of 3.1.2 a, b, c, and d. 

The data- shall be taken 4ising-the same temperatures-and-integration, 
times as were used in the above tets for a total of nine sets of 
data.. 

f. Repeat the~tests^ of 3.1.2 a, b, c, d, and e at three more positions 
of”the ccarse wedge where the positions are set at approximately 
90* intervals. 

g. Re-establish the conditions of 3.1.2(a) and rotate the fine wedge 
at a constant rate equivalent to greater than 0.5s/' sec deviation. 
Continuously- take- data over one complete-wedge rotation. Identify 
the data as: rate tracking data. Data ^hali be reduced and analyzed 
as described in Paragraph 4. 

-3T2“ CARGE-FTEOr POINTING: ACCURACY AND NOISE EQULVAtENT OiSPtACEMENT 

The^object of- ttris test is to determine the relative pointing accuracy over a 
larger portion of the CID. The total range tested is approximately 140 pixels 
in diameter. In addition, the noise equivalent displacement will be measured. 

3.2.1 Setup Requirements 

Rotate the aft and forward wedges until each pair has maximunr deviation as 
indicated by coincidence of the related wedge scribe marks. Identify the pair 
coupled-to the p o t e nt i o met er as fine wedge. 

Establish the conditions of paragraph 3.1.1 b, c, d^ e, and f. 
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3.2.2 Data Requirements 

a. -Obtain tracker position^ wedge ang1e» and intensity data for 60 
positions of the fine wedge, uniformly spaced at intervals of 6* 
wedge rotation. THe tolerance for wedge position settings indicat- 
ed by the pot shall be tO.S degrees and all settings shall be ac- 
quired by rotation in the same direction. Identify the data -as 
■Large Field Pointing*,^ max star,-0*C. 

b. Repeat the test of paragraph 3 .2.2^ (a) above at three more posi- 
tions of the coarse wedge where the positions are set at approxi- 
raatety 30* intervals. 

c. Remove the coarse wedge and repeat the test of 3.2.2(a). 

e. With the CID set at 0*C and the star stationary take at least 50 
position readings at each of the following star intensities: 

Max Star 

BOt 

sot 

32X 

lot 

31 

It 

O.lt 

Jiteatify-the data as-noise equivalent angle- data. 

f. Repeat the- -test of^ paragraph 3.2.2(e) at reduced NORO's as 
follows. Set the integration time for constant intensity readout 
for-each test. - 


3-5 



TR81-04 



36 NDRO's 
16 NDRO's 
4 NDRO'S 
• 1 NDRO'S 

Reduce the data in accordance with the procedure of paragraph 4. 

3-3 ACQUISITION DYNAMIC RANGE 

The purpose of this test is to define the dynamic range for the acquisition 
loop. 

-3-3^U Set-up-and- Test-Requirements- 

With no star input into the tracker, reduce the acquisition and track- thres- 
hold until the tracker acquires on dark current noise, then increase the 
thresholds until it wiTl not remain acquired. 

Increase the star intensity until the tracker refuses to acquire, then de- 
crease it so that it just will acquire. Identify the star intensity as 
"acquisition-max star". 

Now decrease the star intensity until it can not be acquired and identify this 
point -ar^acqusition min-^tar"'. 

3.4 OPTIONAL AND BPLORATORY TESTS 

Additional tests nay be conducted as deemed necessary by the operator. These 
-tests shall not interrupt the sequential taking of a set of data as defined in 
previous paragraphs. 

When such tests are conducted, the purpose and pertinent information shall be 
identified on the data sheets. 


3-6 



TR81-04 



Section 4 
DATA REDUCTION 

This paragraph describes^ the bas^ic reduction processes. Alternate methods- 
that may yield more accurate ro..u1ts may be used at the discretion of the 
personnel involved. 

4.L INTERRIXEL TRANSFER FUNCTION 

A straight line transfer function (outpt vs. displacement) is assuaed over 
each-half-pixel-of-dispTacement. This section defines the constant that must 
be applied to normalize the output to the pixel dimension. 

The ‘four ^•circles*" of 'data taken Tinder paragraph 3.1.2(a) (Max star, G*C) 
shall be used to define the inter-pixel transfer function, 

a. Deviation Wedge Zero Position 

Review each- set-of data, and- tdenttfy. the twa wedge angles at which the 
tracker x-axis output is a minimus and maximum. 

• 

Divide tne difference of these angles by two and identify the result' as 
the x-axis zero wedge position. This is the best estimate of zero devia- 
tion along the x-axis. 

Repeat the process for the y-axis output data and identify the results as 
^-axis-wedge_-Zero. 

b. Transfer Function 


Using- the zero reference position defined in (a) above tabulate x-axis 
outpt vs cos and y-axis output vs cos ^ is the related wedge 

rotat-ion angle form the zero reference. 
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The data cycles through zero,in half pixel Inc r emen ts. Obtain a best fit to a 
straight Tine for each cycle as illustrated below. Identify the integer pixel 
cycle as x and the half pixel cycle as x* as shown.. 


otfTPvr / 



“Extend the lines until they-overlap— two- -zero crossings and determine the 
equivalent differential output for these crossing points. Identify them as ax 
and ax^ respectively for the .appropriate^ cycles- as. shown. Reduce all cycles 
available from each set of data to this form and obtain the mean ax and axf 

The constants for inter-pixel interpretation are determined as follows. 


K = ^ (integer pixel cycle) 
ax 

K , * (half pixel cycle) 

* ax' - - - 

Repeat this process for y axis output data to obtain Ky and Ky'. 
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Multiplying the related output data by these constants wiVl ^ive x and y posi> 
tion in terms of pixels. 

The data obtained In this test shall also be reviewed for systematic non- 
linearities and compared with the previous analytic results that were 
obtained. 

Several -Iterations and alternate- ana^Tytlc methods-may be-used to optimize the 
constants. 


4.2 SMALL AND LARGE FIELH POOfflNS ACCURACY, CIO EFFECTS AND RATE 

TRACKING ACCURACY 

a. Hult1p:lyrtheL fractional .position data of paragraphs 3J. and 3.2 by 
the appropriate constants developed under paragraph 4.1. This 
normalizes fractional position to the pixel dimensions. Perform 
the following analysis on all sets of data taken in paragraph 3.1 
and "372" except for that-of-J.^.2fe). 

b. For each circular set of data above, -*ev1ew the wedge angle at 
which the data was taken and Identify pairs of data points that are 
d1 ametr 1 cal ly“opposed-(approx1mately-180-deg-apart ) _ 

(Jt^l.yi) (xi'yi')..-(x„^ yn^I; CxitA,yni 

c. Compute, for n»l, 2^ 3 n 

' rc«„ - 2 


where R^ Is the deviation of the wedge as measured by the tracker. 
Determine the mean 00 and standard deviation. 
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d. Determine the best center for each circle data points and 
identify it as- (xg, y^). 

Compute: 

for n * 1, 2, 3....n 

Determine the mean 0^,^.) ^d- standard deviation. 

4.3 NOISE EQUIVALENT DISPLACEHENT (NED) 

Determine the mean and standard deviation for the data taken under para- 
graph 3.2.2(e). 

4.4 ADDITIONAL ANALYSIS 


Additional analyses^ and alternate approaches may be accomplished as deemed 
necessary or desirable. However, the analyst is cautioned that alternate ap- 
proaches may result in introduettonjof" 


The greatest contribution to- uncertainty- is- assumed -to— result from- the one 
percent potentiometer. The procedures outlined above are designed to limit 
the resu1tar.t contribution to +0.2S in data uncertainty. 
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-Section 5 

RaATEir EQUATIONS, ANALYSES, AND ERRORS 

This paragraph presents pertinent Information associated with the test process 
and data reduction. 


5.1 DEVIATION WEDGE AND ASSOCIATED-ERRORS 

5.1.1 Deviation Angle 

Figure- 5-1 Illustrates- the_geometry Involved where: 


A » apex angle of the- coirpouond-wedge 
N = Normal to wedge-entrance surface 

1 a the incident angle of ref^ axis (input light axis) relative to N 
r a refractive angle relative to N 
N* a normal to wedge exit surface 
r*" a refractive-arcle-re-lative-to-N' 

1' a angle of exit light axis relatfve^to W 

0 a angle of exit Tight axis relative ^o Input axis (devtatlon-^ngle) 
u a refractive Index of glass 


The baste expression for refraction Is: 
sin 1 

“ *TfiTT 


(5-1) 


and by similarity 

sin 1 ' 
‘ siiT' 


(5-2) 


-f- 

r 




from Figure 5.1 It Is seen that 
r* ^ -(y + A) 

D a 1 + A + r 


(5-3) 

(5-4) 


I 
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By substitution of equation 5-1 through 5-4 we get 

D * 1 + A - sin sin [sin ^ + A]] (5-5) 

The signs for the terns in equations 5-3, 5-4 and S-5 are selected for proper 
polarity for the conditions of Figure 5-1. If the wedge were inverted, such 
as in our rotating.app^lication, the signs of values for (1) and (A) of equa- 
tion 5-5 woald also chanje. 

5.1.2 Error or Uncertainty in Deviation Angle 

From-Flgurfr 5-1 it is seen that the magnitude-of (i) wtlLnot change if the 
wedge is rotated about the reference axis. By similarity, it will not change 
if the wedge is rotated about the normal to the entrance surface (N)^ 

If vfe define the angle between rotation axis and reference axis as (B) the 
angles between the entrance face, (N) and the_rotatTon axis as (o), it is seen 
that for Wie first example [i) would be equal to (o) and for the second 
example (i) would equal (B). 

Therefore, the resultant deviation magnitude (D) would be constant for full 
rotation-of the-wedgei-which— is-the-desire«t-conditton-for^this appTi cation. 

These, however, are unique conditions, where, the two angles did Tiot exist 
simultaneously. 

If the two angles do exist simultaneously, the magnitude of deviation will 
vary as the wedqe is rotated. The^value-of (i)-will tend-to have_a.mean value 
of (o) and will vary from (a + B) to (o - B). The measured deviation (£^), 
will then become 

- - -0- -* D + aO (5-6) 
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Where-(D)-in this case- is the mean devi atiott defTnetTby eg», 5-5, and (AO) the 
error from misalignment. 

It is evident that the effects of (a) and (B) can be included in the overall 
expression for (D^) if their magnitude and direction are known. This however 
would greatly complicate the data handling and evaluation process. 

A simpler approach is to ITmTt the magnitude of (o) and (B) so that the resi^ 
dual error (aD) can be tolerated as-an-uncertainty-in the true devtatlon.- 

The- expression for (aO) ts- foumLby- subst-ltuting^_{o+B) for (i) equation 5-5 
then-ffnding-the difference, in tesu.lts for B=»0 and B^aximum. Therefore, 

Aff = sin"^[p sin-[sTir"^ ^ _ 

(5-7J 

sin'^ [u sin [sin"^ + A]] - B 

which describes the maximum variation from constant deviation as the wedge is 
rotated-. 


S.t."3 Error Budget 


In-our-application-we-are-attenpting-to ver-ffy performance of approximately 
0.4S^ He wllT therefore-arb-itrarily set AO < O.IS. 

The wedge used for the test consists of two 23m wedges each of which is an 
achromatic, crown/flint glass combination having an index of refraction of ap- 
proximately u - 1.53. 

The wedges are rotated relative to each other to provide a fixed apex angle 
that is adjustable from zero to d6n. Since the worse error occurs at maximum 
apex angle we will set A = 46m. 
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The wedges-are - B ou nte d sff^hat the normal to the entrance surface Is typically 
misaligned from the rotating axis by 23n to allow for some assembly tolerance 
we will set a * 30m. 

Substituting these values for &0, u, A and a into equation 5-7 and iterating 
to a solution we get. 



B < 6n 


I for 

AD <-0.1S 

! 

The wedge rotation axis must therefore be aligned to the reference axis, which 
I in this case-wtll be -the star— simulator -axis, to- better than fin to maintain 

the deviation uncertainty within O.IS. 

(T' As the compound wedge an (A) is decreased the tolerance on B increases some- 

what proportionaTTy. 

■ Table 5-1 gives the relationship of B^j^^ vs. A for AO < 0.1s and o < 30a. The 
I approximate mean wedge devition (0}. is also tabulated. 
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Table 5-1 

WEDGE ROTATION AXIS ALIGNMENT REQUIREMENTS (B) 

vs. 

COMPOUND WEDGE APEX ANGLE (A) AND DEVIATION ANGLE (D) 


A 

0 

®max 

46 m 

14G3^s 

6 m 

40 m 

1272 s 

7 m 

30 m 

954 s 

11m 

-20-m 

636 S 

17 m 

10 m 

318 s 

31 m 

5 tn 

159 s 

Slmr 

2itt 

63-s 

:89 m 

1 n 

31s 

13Tm 


Use of the auto collimator feature of the star simulator will result in B< 
2.5m insuring an uncertainty of less than 0.1s. 

5.2 DATA TAKING AND REDUCTIOfT 

5.2.1 Wedge Angle Setting Tolerance 

The test process of paragraph 3 requires taking data at specific positions of 
0 wedc(e rotationv The process results in a. locos Of input^positions that lie 
on a true circle within ±0*1 s (+.0025 pixels) as discussed in paragraph 5.1. 

The radius of the circle (Rn) seen by the tracker, is found by resolving the 
cord length between diametrfcalTy-opposedrdata-points ~as^ described- in para- 
graph 4.2(b) and (c). Since there is an uncertainty in the relative angle 
between data points due ^^pot- readout uncertainty, -a- potent-ial source of 
error exists. 

Assume the data points are taken 180 degrees apart, where is the un- 
certainty in wedge angle between points. The measured half cord length (R^) 
between data points is then; 


5-6 



TR81-04 



180 + 

R„» 0 sin 2 — ^ 

where (0) fs the deviation angle of the wedge. The error resulting, from as- 
soming-the points- are-exactly 180^ degrees apart is therefore: 

, 180 + A*. 

ARn * DLl- sin 2 — *'J 5.3 

If we -wish- to set a limit on A^>to^achieve a given-ARn(max), we get: 

|A*^ai < 2 sin~V- ~^^"^f^V l80 ea. 5.10 

For our application, we will limit the- error to 0.1s, therefore 
ARn(max)=0.1s. Our tests use two deviation- angles. -For- the— small field, 
0<200s and for the large D~1463s. 

Substituting into eq. 5.10, we get: 


Small field A4„,ax < 3.62 degrees 


and 


Large field a^„, 3 x 1-34 degrees 

Since the readout accuracy of the pot is approximately degrees, "the wedge 
rotation positions must be set to accuracies of +2.6 degrees and +0.5 degrees 
respectively for the small and large field pointing tests of 3.1 and 3.2 

5.3 INTENSITY AND DARK CURRENT MEASUREMENTS 

The CIO will be used for all intensity and darK current measurements. If sig- 
nals are converted to equivalent electrons, the dark current JIdj and star 
magnitude are found by: 

T , -4S., , .. 

^d N(lr + 2Ti) sec 
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r 

[ 

i: 

r 





and 


My - 1092^5^. 


4.465 X 10^ N(T^ + 2T^) 
S (f/n)^ 


Where: 


5 dark current e/sec 
S = signal e 

N = Total number of readings In 1st and 2nd set of reads 

Tr = Tine spend taking 1st and 2nd set of reads 

Ti = Integration time between 1st and 2nd set -of reads 

My 2 Star Magnitude 

f/n 5 lens f number setting 


Note that signal from all pixels containing the star must be used for niy 
computation. 


If the digital' s^ignal (Sj) is not converted~to-electrons, using- our scaling of 
179.5 e per LSB we get: 


r 

r 


My - log2ji 


2T48a-xio^N(T;: +-2^7) 


Where Sj Is now the decimel equivalent of the signal directly fron^-the 
tracker. 

5.3.1 Signal Formula 
The basic equation for slyial output Is: 
s - 0.25EsN(Tp ♦ 2T^) 



i 
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-wherer 

s s output signal fn electrons (e) 
s electron generation rate (e/sec) 

N s Total number of NORO's In first and second read 
T,. s total time spent taking the first and second read 
T-f = Integration time between the first and second read 

Figure 5.2 1 nustrates_the_parameters- 



Since the observed signal may be dark current (Ij) or signal current (I5), we 
can use this basic equation for either measurement. By substitution we get: 

T X — ■ 

M n(T + 2TJ sec 
r 1 

and 


^s " n(T + 2T.) - V fi? 

r 1 
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Our process internally subtracts (Tj) when trackings a star, so thatr 


Is " NIT^ + 21^) 

When computing 1^, the signal(s) fr-on alUp-ijcels containing the source must be 
used. 


5.3.2 Star Magnitude Formula 

The signal formula can be used to compute star magnitude by introducing the 
lensirarameters, CID response and star spectral distribution. 

Integration of a GO V, zero magnitude star over GE's published response for 
the CID results in a current of: 

(OMy) Ij = 2.8427 X 10 ® A 

cm 

Where (A) is the combined Itns efficiency and area. Assuming an efficiency of 
0.8, the signal for a given star is: 

. 2 . B427 , 0.8 . 

2.51 '' 


where (r) is the lens radius and. 


■r 


Therefore, 


I 


s 


(f^)^ 

2.51 


and by substitution, 
4S 
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4.465 X 10^ N(T^ + 2T^) 

2.51 S (W 


or 

4.465 X 10® N(T^ + 2T^) Z 
My ■ 1092.51 S" ^f/n^ 


5.3.3 Breadboard FormuTa 

Since our Breadboard lens has an EFL of lO^cnr and-a-variable f/n:_ 

4.465 X 10^ N(T^ + 2T^) 

Our signal is scaled for 179. 5e" per LSB so that if the digital signal (5j,) is 
not converted to-electrons. 


2.488 X 10®N(VtjV 


V ■'092.51- 



[ 

I 


TR81-04 



REPORT APPENDIX B 


FITTING DATA TO THE EQUATION 
OF A CIRCLE BY THE 
METHOD OF LEAST SQUARES 


REPORT NO. SER 86200011 


JUNE 8. 1981 



TRS1-04 



FITTING DATA TO THE EQUATION OF A CIKaE 
BY THE METHOD OF LEAST SQUARES 

A method of fitting data to the equation of a circle by the method -of least 
squares is developed. The method is analogous to a linear least squares fit 
but is valid for circles and can bfr extended witturappropriate caution to :any 
conic section. 

The Problem 


The currently^ planned method for generating well known star position changes 
for testing of the CID tracker system is to rotate glass wedges in the star 
collimator optical path^ This rotation will produce a circular^path-for-the 
focal plane star image. This nonlinear path. presents a problem in data reduc- 
tion. 

The most convient method of fitting the tracker output data to the star posi- 
tion has always been to perform a least square fit of the data to a linear 
scan of the star across the array. A circular star .path requires fitting^ 
tracker data to a circle, and we were unable ta find an-expressioaTor this in 
the literature. • • 

The Approach 

The method of least- squaresrrequires that the sunuof the square of theiTesi- 
•duals, vi» in the- following expression be minimized. 

V- * y. - y. 

where y^- is the value of y obtained by substituting x-j »- x into the- equation 
y - f(x). 

If— the-equaticn is-of the form: 


1 
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f (x. y) - 0 


Vi * f (xi. yi) 


If we define this sum as S we may say: 


S 5^ 


i 


n 

I V 
« 1 


2 

1 


The expression, f (x), will have-unknowns; aj -ajj ^ich are to be evaluated^ 
witt the n, (xi, yf) sets. The necessary and sufficient condition which 
minimizes S is that: 

„ fiS . «S ^ 

71“ = 7a" * 0* 7a ■ 0* 

1 ^ k 

A circle can be expressed as a polynominal in two variables, f (x, y) » 0. 

(x - xq>^ + (y - ^ • 0 

or 

x2 + y2 - rxg x - 2 yo y + xo^-*-yo^ — “ 0 


-Me know that the only unknowns are xO *"<1 ^ 0 * location of the center -of 
the circle and R:, its radius. We will set: 


ao = -2x0 al » -2yo 


and 

*2 ■ xo^ 

It is at this point that our analysis ceases to be entirely rigoris. The 
partial differential equations: 
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I 
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1 

I 

r 




are required to be independant of one another, and it is dear that in this 
construction a2 is a function of ai and ag* 


We can extract ourselves from this dilema if we have an approximate knowledge 
of the location of the center of the circle. If this approximate value is 
subtracted from the data such that ^2 • independance of the partial 
differential express^ions is preserved and we may-proceed_wlth our developments 


Using the substution above our expression for a circle becomes: 




Remembering that: 

Vi « f(xi. yi) 

we have: 

n - n , 

S 2 r V.2 » r f(x., y 
i » 1 ^ i » 1 ' ’ 


and 

I 


? 1 2 ‘i * “O "l * "l * *2> ■ 

% ? 1 2 y,(x,2 . ,,2 • a, X, * y, a a^) - 0 

■ %* J 1 * *2’ ■ “ 


3 



TR81-M 




1 , 

!'. 

1 , 

i 

I 

I 

I 

I 

I 

1 

\ 



abbreviating S to r. and dividing all three equations. by two we obtains 
i » 1 

r + t Xi yi^ + ao t Xi^ + ai_t xi yi + az x xt - 0 
X xi^ yi + s yi^ + aO x xi yi + al x yi^ + az x y^ ■ 0 
E X.2 + X y^2 + j. j + az n - 0 

In matrix form this may be expressed as: 


— 




r 3 , 3 

“o 


Ex/ + X X. y^ + X x^ 

- 

-xx^ -X x^y^ 


x 

c X, y, a I y,^ * I yi 

m 

-t x,2 y, -.1 y,3 

J2 


Ex + X y. + n 


-IXi^-tyr/ 


After solution of this expression for ao. al, and az^ 


and 



*0 *1 , 

- “4 + -y- - >Z 


The residuals are evalulated by substituting the input data pairs into the 
original oxpression aftar the a» l>a»e >*»» a«alnate<l. 

■ xi^ ♦ yi ao xi ♦ al yi + a2 
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As in a linear least squares analysis, the standard deviation of the v. the 
nx}st reasonable measure of the quality of the f tt^ 

Applicability 


It should be noted that this method is useful for any two dimensional expres- 
sion in which the power of the variables take on positfvr integer values -ds 
long as the data can be manipulated to maintain independence of the partial 
differential equations. 

The general expression for a conic section is clearly applicable: 
bo + bi y^ + b 2 X + b3 y + b 4 “ 0 
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